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We introduce a group theoretically motivated procedure of parametrizing non-forward matrix 
elements of non-local QCD operators by (two- variable) distribution amplitudes of well-defined geo- 
metric twist being multiplied by kinematical factors (related to the Lorentz structure of the operators 
and to the target states) as well as position-dependent coefficient functions resulting from the (in- 
finite) twist decomposition of the operators. These distribution amplitudes are interpreted as (sum 
over) power corrections of the double distributions. — Using the technique of harmonic polynomials 
for the local operators we determine the (infinite) twist decomposition of totally symmetric opera- 
tors completely and for operators with non-trivial symmetry type up to twist r = 3. This covers 
the phenomenological interesting quark-antiquark operators. Using these results we determine the 
power corrections to the various double distributions and the vector meson wave functions. It is 
shown that the structure of the kinematical power corrections may be obtained, by harmonic exten- 
sion, from the corresponding expressions for operators or distribution amplitudes, on the light-cone. 



PACS number(s): 12.38 Bx, 13.85 Fb 



I. INTRODUCTION 



The universal, non-perturbative parton distribution amplitudes parametrizing, modulo kinematical factors, the 
matrix elements of appropriate non-local quark-antiquark (as well as gluon) operators play a central role in phe- 
nomenological considerations. These operators occur in the quantum field theoretic description of light-cone dom- 
inated hadronic processes via the non-local light-cone expansion For deep inelastic lepton-hadron scattering 
and Drell-Yan processes the parton distributions are given as forward matrix elements of bilocal light-ray operators 
resulting from the time-ordered product of appropriate hadronic currents. For (deeply) virtual Compton scattering 
and hadron wave functions the so-called double distributions and hadron distribution amplitudes, respectively, are 
given by corresponding non-forward matrix elements. Thereby, various processes are governed by one and the same 
(set of) light-ray operators. 

For the sake of definiteness let us remember the amplitude of virtual Compton scattering, 

T^AP^,Q■,Si) = i J d^xe"'^P2,S2\RT{J^{x/2)M~x/2)S)\Pi,Si), (1.1) 

where = —q^, q — <}2 — Qi denotes the momentum transfer and S is the (renormalized) S"— matrix. Restricting to 
lowest order the renormalized time-ordered product of the electromagnetic (hadronic) currents, Jy.{x) =:'0(a;)7^?/'(a:) : 
in the vicinity of the light-cone, x ^ x, Sp' ^ 0, where 

X ^ X - ciixQ) - {xCf - x"^) with C^ = l, (1.2) 

is approximated by (a sum over) coefficient functions cju)a{x'^ , kx, fi^) times the matrix elements of (renormalized) 
light-ray operators (/ denote the quark flavours) 

(P2,^2|Oi5^(«i,-«i)|Pi,^i) = V^(P2,^2|V5^(0)(75)7a(iB)"^/(0)|Pi,^i). (1.3) 

n=0 
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Similar considerations hold for the (vector) meson wave functions occurring in processes like exclusive semileptonic 
(or radiative) B decays and hard electroproduction of vector mesons {V = p, w, K* , (p). The physical interest in this 
case arises from the fact that the following vertex function is directly observable in QCD: 



I 



d^xe''i^0\RT{J^{x/2)J,{~x/2)S)\V{P,\)), 



(1.4) 



where can be axial or vector currents and \V{P,X)) is the vector meson state of momentum P and helicity A. 
Here, the universal nonperturbative quantities, the meson distribution amplitudes (DA), sometimes also called meson 
wave functions, are given by means of the vacuum-to- meson matrix elements of bi-local operators which, in the limit 
X ^ X, read 



with r = {l,75;7a,7a75;0-a/3,a-a/375}- 

Possibly, in the near future the experimental precision data will allow for the determination of non-leading contri- 
butions and, therefore, require for a careful analysis of the various sub-dominant effects contributing to the physical 
processes. When considered beyond leading order, i.e., beyond lowest twist operators in tree approximation, one is 
confronted not only with radiative corrections but also with power corrections resulting from higher twist as well as 
target mass effects. Higher twist contributions are obtained by the decomposition of the bilocal operators Ot{kx, —kx) 
with respect to (irreducible) tensor representations of the Lorentz group having definite twist r = (scale) dimension 
minus (Lorentz) spin . These representations are characterized by their symmetry type (under index permutations 
of the traceless tensors) which is determined by corresponding Young frames having, in the case of the Lorentz group, 
at the most three lines, [to] = [TOi,m2,TO3]. 

The twist decomposition of non-local light-ray operators ('finite on-cone decomposition'), as far as they are relevant 
for hadronic processes, has been completed quite recently. A comprehensive presentation of these results can be found 
in Refs. where also earlier studies are mentioned. Recently it has been used for the definition of quark distributions 
functions of well-defined geometric twist together with their unique relationship to the quark distributions of dynamical 
twist (being introduced in Ref. Q). However, in order to get information about the target mass contributions one 
is forced to consider the twist decomposition off-cone, taking into account all the trace terms leading to expressions 
suppressed by powers of M^/Q^. 

For local operators which are characterized by symmetry type [n], i.e., being determined by totally symmetric 
traceless tensors, there exists a well-defined group theoretical framework |^ which works on the light-cone as well 
as, by unique harmonic extension, also off the light-cone ('infinite off-cone decomposition'). For any other symmetry 
types which, however, are relevant in the above mentioned physical applications no general theory exists. But, from 
our earlier papers partial results can be obtained. Using them we are able to present the decomposition of the 
relevant off-cone operators up to twist t — 3 which already covers the cases being experimentally accessible in the 
future. In principle, twist-4 results also would be of interest, but there occur at least two bi-local operators of different 
symmetry types of which only one has been fully determined. 

Let us now explicitly list down, suppressing their fiavour content and, as usual, the gauge link, the non-local quark 
operators which will be studied in detail. We determine the twist decomposition of the following chiral-even (axial) 
vector operators 



(0|Or(i,-i)|F(P,A)) = ^(0|V5^(0)r(xi^)V'(0)|V^^^'^H^,A)), 



(1.5) 



n=0 



Oa{KX, —Kx) = ■li^{Kx)ja1Pi — Kx), 
Oza{liX,-Kx) = ■lp{Kx)jal5'>Pi-Kx), 

together with the corresponding (pseudo) scalar operators 

0(5) (kx, — Ka;) — x"0(^5ja{Kx, ~Kx), 
and the chiral-odd scalar and skew tensor operator 

N{kX, —kx) — 'lp{Kx)'lp{—Kx), 
M[af3]{KX,-Kx) = ■4)(Kx)aap^{-Kx). 

together with the vector and scalar operators 

Mainx, —kx) ~ x^ M[g^/j^{Kx, —kx), 
M{kx, —kx) = x''9"M[q^] (kx, —kx). 



(1.9) 
(1.10) 



(1.11) 
(1.12) 



(1.6) 
(1.7) 



(1.8) 
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The operators ( |l.6| ) - (1.12) constitute a basis not only for the (usual) parton distributions as well as meson wave 
functions [p|~|ll| but also for the consideration of double distribution amplitudes being relevant for the various light- 
cone dominated QCD processes under consideration. Here, their twist decomposition off the light-cone will be given 
up to twist 3 and, in the case of scalar operators, also for any twist. In fact, the 'external' operation of contracting with 
or x^d"" also influences the possible symmetry type of these operators and, therefore, of their twist decomposition. 
Notice that the external coordinates are not multiplied by k. 

As a result these off-cone operators of definite twist are given for the local as well as the resu mme d nonl oca l operators. 
The local operators given as the n-th Taylor coefficients of the non-local ones, see, Eqs. (O) and ( |l.5| ) (or, more 
generally, Eqs. (2.1) and (2.2) below), are represented by (a finite series of) Gegenbauer polynomials C^(2:), u > 1. 
The nonlocal operators, being obtained by resummation with respect to n, are represented by (a related series of) 
Bessel functions or, more exactly, either J^_i{z) or I^_i{z) depending on the values of their arguments. 

The group theoretical method for the determination of target mass corrections in unpolarized deep inelastic scat- 
tering using harmonic scalar operators of definite spin and the corresponding matrix elements in terms of Gegenbauer 
polynomials has been used for the first time by Nachtmann |]l2| . Some years later, using the same procedure, the 
target mass contributions for polarized deep inelastic scattering were studied in Refs. ||l|-|l|l. A short review of 
Nachtmann's method in deep inelastic lepton-hadron scattering was given in Ref. [0. This method differs from the 
one being applied by Georgi and Politzer ||l^ for the target mass corrections in unpolarized deep inelastic scattering. 
The latter method which is tailored to the forward case has been applied also to the study of target mass corrections 
of polarized structure functions . 

Here, we generalize Nachtmann's procedure to the case of non-forward matrix elements applying it to off-cone quark- 
antiquark operators. In this paper we restrict ourselves to a fairly complete consideration of the various distribution 
amplitudes and postpone a more general application, e.g., to the virtual Compton scattering (using a somewhat 
different approach the latter has been already considered in Ref. |^ll). However, we determine completely the off-cone 
power corrections of the meson distribution amplitudes in a;— space which are much easier to handle. 

The paper is organized as follows. In Sect. II we present our method by considering a typical example. The power 
of that approach consists in the fact that one determines the twist decomposition at first for the local and nonlocal 
operators and only afterwards takes the matrix elements. In Sect. HI totally symmetric operators are studied thus 
obtaining the mos t gen eral result including the whole series of infinite twist. In Sect. IV the twist decomposition of 
operators (1.6) - ( 1.12 ) up to twist 2 and 3 is given. In Sect. V these results are applied to the double distributions 
and the meson distributions. 



II. NON-FORWARD MATRIX ELEMENTS OF NONLOCAL OFF-CONE OPERATORS: THE METHOD 
A. Parametrization of non-forward matrix elements by independent double distributions of definite twist 



To begin with we consider the bilocal off-cone quark-antiquark operators (1^), ( |l . 7| ) , ( |l.9| ) and ( 1.10 ), i.e., operators 
without external operations which, generically, will be denoted by Or{Kx, —kx). More generally, they are given as 

Or{y + Kix, y + K2x) = ^{y + nix)TU{y + kix, y + K2x)ilj{y + K2x), 

with rj = y + {k2 + ki){x/2) and ^ — 2kx = (k2 — ki)x being the centre and the relative coordinate, respectively. 



U{y + Kix, y + K2x) — V exp | — ig 



being the (straight) path ordered phase factor ensuring gauge invariance (which is omitted in the following). 

The Fourier transforms of the 'centred' operators Or{nx,—K,x) and their nth moments C'rr!(a;), i.e., their Taylor 
coefficients w.r.t. k, are related as follows: 



Or{Kx,-Kx) = / d\C'r(<?)e'''^« = ^ 



(i«)^ 



(2.1) 



n=0 



with 



Ovn{x) = / A\Or{q){xqr = {-\T^OT{nx,-Kx) 



(2.2) 
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For notational simplicity we wrote the Fourier measure without the usual factor l/(27r)^; obviously, q should not be 
confused with some momentum transfer. It should be remarked that the Taylor expansion of the non-local operators 
into local ones is only justified in a restricted region of the Hilbert space Here, we assume their existence after 

taking physical matrix elements. 

Now, we introduce a suitable, but preliminary parametrization of the non-forward matrix elements of the bilocal 
operators Orinx, —kx). Namely, they can be represented as (everywhere, we use summation convention w.r.t. a) 

(P2,^2|Or(«:x,-«x)|Pi,5i) =/C?(P) j DZ e'«(^^)^ /„(Z, P,P„ x^; (2.3) 

This parametrization which takes up and generalizes previous ones P, |2^ , p3[ | deserves some explanations and further 
comments: 

• First, we introduced the notation P = {P+,P_} and Z = with P± = P2 ± Pi and z± = \{z2 ± zi) 
thereby defining some (2-dimensional) vector space with scalar product PZ = ^P^z^ = P+2;+ + P-z^. In 
addition, the integration measure is defined by DZ = dzidz20{l — zi)6{zi + 1)9{1 — Z2)9{z2 + 1). 

• Next, /Cp(P) denote the linear independent spin structures being defined by the help of the (free) hadron wave 
functions and governed by the F— structure of the corresponding nonlocal operator Or- For example, in the case 
of the virtual Compton scattering, there are two independent spin structures, the Dirac and the Pauli structure, 
/Ci = S(P2, S2)lMPi^ Si) and ICl = u{P2, S2)a^,,P1u{Pi, Si)/M , respectively. 

• Furthermore, modulo these spin structures, the Fourier transforms /a(Z, P^Pj, a;^; /x^) of the matrix element 
(P2, S2\Or{Kx, —Kx)\Pi, Si) with respect to the independent variables k{xP) are the associated (universal, renor- 
malized) two- variable distribution amplitudes as introduced in |^ but here extended off the light-cone. As it is 
obvious, also their dependence on PiP,, x'^ and the rcnormalization point /i^ has to be taken into account. 



The representation ( ^.3| ) gives the most general expression for the non-forward matrix elements of Or{KX, —kx). 
An essential aspect of that approach consists in the support restriction, — 1 < 2^ < 1, of the distribution 
amplitudes /^(Z, P^Pj, a:^; /i^) since the matrix elements can be shown to be entire analytic functions with 
respect to xPi (for a detailed discussion, see, 0,0|). 



Obviously, the representation ( p.3[ ) may be used also for more involved matrix elements, like (scalar) meson 
production, by extending the final state to (Pj , 6*2 ; fci | thereby enlarging the Z— and the P-space and extending 
the set of possible spin structures K^(F) and related distribution amplitudes fa (see, e.g., Ref. Q). It can be 
specified also to vacuum-to-hadron transition amplitudes. 

Finally, let us comment on operators with external operations, i.e. when the (axial) vectors or the (skew) tensors 
are multiplied 'externally' by x^ and/or d". In these cases the 'external' vector x^ or tensor x^ d" i s assumed not 



to be Fourier transformed and the external operations have to be applied onto both sides of Eq. (2^). Thereby the 
tensor structure of these external operations matches with the tensor structure of which by itself are independent 
of the coordinates. Despite of this the external operations heavily influence the possible symmetry type of the local 
operators and their decomposition into irreducible tensor representations of the Lorentz group. 

Now, let us take into account that the non-local operators Or^nx, —kx), formally, are given by infinite series of 
operators of growing (geometric) twist r (using summation convention also w.r.to F), 

Or{Kx,-Kx)= c^''f{x)0^fMKx,~-Kx), (2.4) 

T>T„i„ 

with 

o'f \Kx, -Kx) = v'-^'^^'Or'iKx, -Kx), (2.5) 
(pir) X V^^">)/' = S^^'V^P^', (2.6) 
where Ppp^/ are well-defined projection operators. Also here we have to add some comments in order to make the 



content of Eqs. ( p.4[ ) - ( p.6[ ) more definite 

• The projection operators V^^^^ {x, dx) which immediately act on the undecomposed operators Or{KX,—Kx) 
depend on the coordinates and their derivatives and, eventually, contain additional integrations with respect to 
some auxiliary variables. Usually, the summation over F' is restricted to the same tensorial type. The specific 
form of the projections may be read off from the explicit twist decom p ositio ns, cf., e.g.. Sec. Ill, where various 
infinite twist series have been derived, as well as Sec. IV, Eqs. ( 4.17 ) - ( 4.19| ). 
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The coefRcient functions c"'^^ (x) essentially depend on (powers of) x'^ but, in principle, according to the tensor 
structure of the operators they also could depend on x'^, cf., e.g., Eqs. ( [3.19 ). In addition, we remark that 
different operators of the same twist may occur having, of course, different coefficient functions. 



• Approaching the light-cone, a; — > i, the series ( p. 4] ) terminates at some finite value Tmax since almost all of the 
coefficient functions vanish for = 0. That situation has been considered extensively, e.g., in Refs. 



O bvio usly, analogous decompositions exist for the local operators Ornix) which are uniquely related to Eqs. (2.4) 
- (2.6). They have been derived in parallel to the non-local operators in Sees. Ill and IV. 



• Finally, we remark that the relations (2.1) and ( p.2[ ) between the local and the non-local operators, Ornix) and 
Or{KX, —Kx), and the Fourier transforms, Or{q), also hold for the corresponding operators of definite twist r. 

Introducing the decomposition ( |2.4| ) into the representation ( |2.3| ) allows for an analogous decomposition of the distri- 

\F,,x^^), 



E 

'^^'^mi 



c(^r(2^)/C?,(P)/M(Z,P,P„x2;/i2). 



(2.7) 



bution amplitudes, /a(Z, 

/Cf(P)/a(Z,P,P„a;2;^2) ^ 

(r) 

Then, for the matrix elements of operators Op with definite geometric twist r we obtain two different versions, 
{P2, S2\0'f\Kx,-Kx)\Pi, Si) =/C?(P) J DZ e^'f^'^^^/f '(Z,IP»IPj,a;2;^2-) 

= r^^f{x,d,)IC^,{F) f DZe-(-')Vi"n^;A''), (2.8) 



where the first line is only the red uctio n of relation (2.3) to their twist components and where, in the second line, using 
properties ( |2.6| ) of the projection (2.5) onto the non-local operators of definite twist we introduced the related Lorentz 
invariant double distributions fa'^\'Z; /i^) of definite twist t not suffering from any power corrections. Concerning that 
definition of double distributions let us give a heuristic justification and some of its consequences: 



• Reminding the Fourier representation (2.1) for the non-local operators Oy\kx, —kx) it looks reasonable that the 
x^— and P— dependence of the distribution amplitudes /a(Z, PiPj, x^; yit^) should be already uniquely determined 
by the projection operators of ( |2.5| ) acting on the exponential on the right hand side of (2.8) and, possibly, on 
some external operations matching /Cp, (P). This reflects the fact that every local operator of definite twist 
is determined by its symmetry type expressed, within the polynomial method, through specifi c diff erential 
operators acting on (traceless) harmonic tensor operators containing the a;^— dependence, cf., Eq. ( ^.ICj ) below. 



As a consequence, the mass corrections of the matrix elements (2.3) are determined b y E q. (2.8). Likewise, 
the distribution amplitudes /a(Z,PjPj,a;2;/i2), using both forms of the matrix element ( |2.8| ), formally may be 

expressed through the double distributions /^^■'(Z;/!^) according to: 

/M(Z,P.P.,^';A*') - {lC-\¥))l(e~^^(^^^^V^p^\x,d.)e^^^^A IC^,{P)fi:\Z;fi') 



^ Ti:}{{x¥Z),x'{FZ)') fir\z■,^^^ 



(2.9) 



here, jr^^?((xPZ), (PZ)^^;^) contains any information about the power corrections of the double distribution. In 
addition, we observe that, as a result of the twist projection, the (xP)— dependence of the matrix element is 
only partly determined by the double distributions fj^j\Z; ^'^). 

• Obviously, since the decomposition of the distribution amplitudes is uniquely related to the decomposition of 
the corresponding (quark-antiquark) operators the power corrections of the double distributions, Eq. (|2.9D, are 
already determined by the twist decomposition o/ these operators. Both decompositions, besides on (xPZ), also 
depend on the following combination of variables, (xPZ) j \J x^ (PZ)^ and ^ (xPZ)^ — ^^(PZ)^, being the argu- 
ments of appropriate Gegenbauer polynomials and Bessel functions in the case of local and non-local operators, 
respectively; cf., Subsecs. II B and II C below. 

• Notice, that the double distributions of definite twist fi'^\Z; fj?) being independent on x"^ are already determined 
by the decomposition of the non-local matrix elements on the light-cone. Also their renormalization properties 
are determined by the light-cone operators only. 
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This procedure of introducing double distributions of definite geometric twist generalizes to arbitrary off-cone values 
of a; what has been introduced in for the parton distribution functions of definite geometric twist, i.e., for the case of 
forward matrix elements restricted to the light-cone, and later on applied to the (vector) meson distribution amplitudes 
[01, i.e., to vacuum-to- meson matrix elements, also restricted to the light-cone. 

Let us now present the arguments leading to the r.h.s. of Eg. ( ^.^ ) more explicitly. Namely, the non-forward matrix 
elements after performing the twist projection of relations {2A) are given as follows: 

{P2,S2\0\:\Kx,-Kx)\P,,Si) = {P2,S2\V^^f{x,d) J d4gOr'(9)e''=^«|Pi,5i) 

= ld\(v^^f{x,d)e"^^'^){P2,S2\Or'{q)\Pi,S^) (2.10) 

Dzh^^f'ix, d,) e-(-P)z) JC-, (P) /(-)(Z; m'). 



Here, in the second line we simply reordered the operations of integration, taking matrix elements and twist projections 
appropriately. In order to get the third line we observe the symmetry in x and q of the local twist projections, 
'^nr' d){xq)" = V^^^Y dq){xq)" , from which one derives for the second line: 

E ^ {^nv\^^ d) X^^... X^-) I d\ {Vi^fiq, d,) q,,... q,„){P2,S2\Or-'niq)\Pl, S,). 



n=0 



Here, the integrand contains the matrix elements of the irreducible local operators of definite twist o'f^^ ^ (q). Their 

reduced matrix elements fri\n2, J^i + ^2 — ^1 which are related to the decomposition of {P2 + Pi)" into independent 

monomials may be represented by (double) moments of corresponding double distributions f!l'\zi, Z2) multiplied 

by the kinematical factors /Cp(Pi,P2)- Finally, since q reflects the dependence of the (local) operators on D, after 
performing the integration and the resummation over n, q within the integrand simply has to be replaced by PZ. 
That procedure will be demonstrated explicitly in Subsec. H.C. 

Formally, in the various experimental situations, in any Fourier integrand having the structure of Eq. ( 2.1C ), we 
only have to perform the following replacements (denoted by =): 

(A) In the case of non-forward scattering, as just explained, we obtain 



(F2,52|Or(g)|Pi,5i) =/Cf(P) J DZ<5(4)(g-PZ)/M(Z;A*'). (2.11) 
(B) In the case of forward scattering, i.e., for Pi = P2 = P, the situation changes into 

{P,S\OriqmS)=ICUP) I dz5'-'\q-2Pz)f^j\z-^?), (2.12) 



with _P+ — 2P, P_ = 0; z+ = z. These distributions f^\z,^'^) are obtained from the double distributions f^\'L, fj?) 
by integrating out the independent variable z^, i.e.. 



/M(z,m') - j dz_/M(z+ =z,z_;/.2). (2.13) 
(C) In the case of vacuum-to-hadron transition amplitudes, e.g., for the meson distribution amplitudes one obtains 

{0\Or{q)\V{P,X))^jCUP) [ d^S^'Hq-POfY\^;t^'), (2-14) 



with P2 = P; -^2 = C- Obviously, the spin structures are different from the above cases and, by construction, only P 
and ^ occur. 

From the foregoing discussion it is obvious that the mass corrections to the various physical processes are completely 
determined by the twist structure of the bi local op erato rs Oy\kx,—kx) = J d'^qO^\q) e^'^^'^ . Performing matrix 
elements according to the replacements Eqs. (2.11) - ( [2.14 ) we obtain the related expressions for the 5-matrix elements 
of the physical process under consideration. Then, a further Fourier transformation with respect to the coordinate x 
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like in Eqs. (1.1) and ( |l.4| ), finally, leads to a representation of the scattering or transition amplitudes in terms of the 
inverse momentum transfer depending, quite generally, on PiPj/Q^ . 

In the next two subsections we demonstrate our procedure for the simplest nontrivial case, first, directly for non- 
forward matrix elements and then, more efficiently, for the corresponding nonlocal operator itself. Thereby, we 
introduce the technique of expansion into (local) harmonic polynomials and their resummation to non-local harmonic 
operators. 



B. Non-forward matrix elements of bilocal operators having definite twist: A simple nontrivial example 

Let us now demonstrate our procedure by the simplest nontrivial example, the non-forward matrix element of the 
twist-2 part of the (pseudo)scalar operator, Eq. (1.8), which is relevant for the leading terms of virtual Compton 
scattering [p5| . This operator containing the 'external' truncation by is given by [p6|,p| 



0''"^{Kx,-nx) = i^{Kx){x-f)^{~Kx)+Y.J^ dt(^)' '(^) \,(°_ 1), ^(^fa)(a:7)V'(-^to). (2.15) 

1 matrix elements, taking into account the Dirac structure only, read 
{P2\0''"^{KX,-Kx)\Pi)^u{P2){x-f)u{Pi) [ DZe"=("^)^F^'^(Z) 



oo p\ ^ f k—1 2/^1 — \k p 

+ i:J^ dt(i^) " (^) ^f(^ Hm^lHPi) J DZe-(-)-Fi^)(Z), (2.16) 

where F^\'L) denotes the twist-2 Dirac type double distribution. Obviously, on the light-cone only the first term 
survives. Here, we take the complete expansion into account. Performing the differentiations. 



IKtxVZ 
^ 1 



□fc(2,^)ei«:tePZ ^ i2k (itixmf^ + 2fct2''-i (iK7PZ)(ma;PZ)2('^-i) 

and using 

(1 - tf-H'+^ ^ ^ (i«xPZ)" (n + £+l)! 

(fc-1)!^! ^ n\ (n + k + £+l)W: 

^ ' n — 

this can be rewritten as follows: 



n\ (n-f 2fc)!(fc- 1)! 

Now, shifting fc fc + 1 in the second fc— summation, and observing the series representation of the modified Bessel 
functions of the first kind (cf., Ref. H, Eq. 1.5.2.13.27), 

E ^"^^"^tn ■ ^" - ^"(^Z^) (2.17) 

^-^ 1 n + 2i^ + 1 n! 

n— 

for V ^ k + 3/2 and z — iK{x¥)'L we arrive at 

{P2\0''"'^{nx,-nx)\Pi) = j DZi^^^^(Z){u(F2)(.T7)"(^'i)(2 + xa)- |ira2M(P2)(7ra)"(A)}(3 + a;a) 



X 



^ {{Kxf{nf/Ay ^ (n + fc+1)! (iKa;PZ)" 

k=0 n=Q ^ ' 



DZ F^^^ (Z) u{P2)-^^'u{Px)\x^, (2 + xd) -\\k Y^TL | (3 + xd) 

E ((^^)'(^^)V4)' (i^^pz)-fc-3/2 4^3/2 (i^^pz/2). 



fc! 

fc=0 
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Now, using I^(&''/'^z) = J„{z) with the (usual) Bessel functions J^(z) as well as (of. Ref. [pTl, Eq. II.5. 7.6.1) 



fc=0 
„2^2 , 



for 



|2t| < z, 



(2.18) 



and substituting t = kx^ / {Axq) and z = Kxq/2 we finally arrive at 

{P2\0^'^^{nx, -Kx)\Pi) =^u{P2)-i^u{Pi) fBZFj^^Z) {.T^ (2 + a;9) - |iKP^Za;2|(3 + xO) 



(kV(2;PZ)2 -x2(PZ)2) ^^^3/2(1 \/0rfraP~^2?2(PZ)2) 



(2.19) 



This ho lds in the case (xPZ)^ — a;2pz^ > 0, otherwise we have to change into Iv{z)- That result exactly corresponds 
to Eq. (2_^) for the example under consideration. Observe that another half of the exponential is concealed within 
the Bessel function! 

In the case of forward scattering because of u{P)j'^u{P) = we obtain 



{P\0''"^iiix,~Kx)\P) =2y/^ I dzF''^\z)\^{xP){2 + xd) 



IKZ P^X^ 



Ys + xd) 



(2Kz^/{xpy 



-3/2 



J- 



3/2 



(^Kz^/{xpy 



(2.20) 



If this expression is restricted to the light-cone, 0, the well-known parton distribution is obtained. Namely, 

using the Poisson integral for the Bessel functions (cf., Ref. |28|, Eq. II. 7. 12.7), 



r(.-^i)J.(z) = ^(J) / dt(i-iY-^/^e"^ 



/^V2 



for Re > — 2, 



(2.21) 



we obtain (^ = kz{xP)) 

{P\0''"'iKi,-Ki)\P)^2{iP)V^ J dzFl^\z) {2 + Cd^)i3 + Cd^){2(y^^'jy^{Oe^^ 



2{S:P) I dzF^^\z) I dtt{l-t){3 + tdt){2 + tdt)e^'''^ 



= 2(iP) J dzF^j^\z) dti(l-Htat)e2'*« = 2(iP) j dzF'j^\z)e 



2iKz(£P) . 



(2.22) 



in the second line we introduced the Poisson integral after shifting the integration variable t ^ (t -|- l)/2 and followed 
by changing — > tdt^ then we partially integrated two times retaining finally only the surface term at t = 1. As a 

result we arrived at the twist-2 parton distribution as introduced in Fj^\z) = F^'^''{z), which coincides with fi{z) 
in the notation adopted by Jaffe and Ji . 

The power corrections of the parton distribution F(2)(z) are obtained according to the relation (O). Namely, using 
the matrix clement ( ^.20 ) and, again, taking into account the representation (2.21), one gets {Y = -^/l — x-^P^ /{xP)"^): 



F'^^>{z,xP,x^P- 



■ ^2) - F^^\z) e-2-(-^)-l|(2 -f xd) - iKz f-^|(3 + xd) f\t (1 - t^) e''^* 

8 1^ (xP) J J_i 

= F(^)(z)e"^'''(^^)^ J_|(i + y)2giK(^P)^(i+>') _ (1 - y)2e'''(^-^)''(^~'*')| 



(2.23) 



Analogous expressions result for the pseudo scalar case by observing u{P)j^j5u{P) — 2S^. 

The situation in the non-forward case is somewhat more complicated and, therefore, will be omitted here. There, 
not only th e Dir ac but also the Pauli structure has to be taken into account whose general form coincides with the 
expression (2.19) wit h /C^ being replaced by /Cp but, as can be proven by using the hadrons equation of motion, in 
accordance with Eq. (2.£) both the Dirac and Pauli structures mix with each other. 

Finally, let us point to the fact that expressions similar to Eq. ( ^.2C| ) have been obtained in Ref. |2^ for inclusive 
particle production in e+e^-annihilation. This work is based on a procedure | p6[ | which, in order to obtain higher 
twists, avoids the consideration of group representations and, instead, makes use of the quark (and gluon) equations 
of motion. As already mentioned |^ in the case of vector (and tensor) operators by this approach the tracelessness of 
the (local) operators and, consequently, also the definiteness of their geometric twist may be missing. 
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C. Nonlocal off-cone operators of definite twist: The technique of harmonic polynomials 



Obviously, the a;— dependence of the expressions ( ^.19 ) - ( 2.22| ) is completely determined by the twist structure 
of the operator from which the matrix element has to be built. Therefore, let u s der ive these results directly from 
the related twist-2 operator. But, instead starting from the nonlocal expression ( 2.1^ ) we use its local version being 
given by harmonic polynomials of order n in x. This polynomial technique uses the vector x S M'* as a device for 
writing tensors with special symmetries in analytic form [8 pO|j3l[ |. It has the advantage to be directly related to the 
irreducible tensor representations of the Lorentz group. Its group theoretical background as far as it is related to 
totally symmetric tensors has been given in Ref. B . 

The local scalar operator according to (2.2) is given by 



0,,+,{x) = 4>mx-f){xDri;{0) = / d\{{pri^){q)x^{xqr, 



(2.24) 



which shows that, in this connection, q formally replaces the covariant derivative D sandwiched between the quark 

operators^ This local operator has a (finite) twist decomposition whose complete series will be given later on, 

cf., Eq. (3.11). Its twist-2 part is simply given by the harmonic polynomials l32,Hl, 



[^] 

E 



{-l)''{n+l - k)\ 



4'=fc!i 



n 



D! 



n+l 



{x)^vl^l{x\d^)0^^,{x), 



(2.25) 



being characterized by traceless, totally symmetric tensors of rank n+l whose indices are completely contracted by 
x^^i . . . a;''"+i. They obey the condition of harmonicity: \30^^_^i{x) — 0. 
Performing the derivatives, 



x"'n''{-fx){qxY 



n 



i-fx){q'x'nqx) 



\n—2k/ 



2k nl 



{n-2k)r ' '^■"^ """^'^ ' (n+l-2fc)! 
and using the series expansion of the Gegenbauer polynomials (see, e.g., Ref. [p7|. Appendix 11.11), 



x\^q){q'x'f-\qx) 



-l-2fe, 



(2.26) 



{-l )''{n-k + u- 1)! 
klin-2k)\ 



i2zr 



-2k 



we obtain 



1 



d\(^7'^^)(g)L^(|y 



q'^x^ 



C„ 



qx 



c,. 



qx 



vV 



(2.27) 



(2.28) 



Again, some remarks are in order: 



• The factor l/(n+l) in front of the integral results from the equality a;p(a;g)" = {l/n+l)df^{xq)'^~^^; it corresponds 
to the normalization of the decomposition of the local operators into irreducible representations, i.e., to the 
nor maliza tion of the associated Young operators. In the nonlocal formulation it leads to the i— integration in 
Eq. ( ^l6| ), cf., also, Ref. §. 



• The (discrete) 6*— functions guarantee that both the fc— summations do not contain undefined factorials thus 
correctly leading to the corresponding Gegenbauer polynomials. In the following this is circumvented by the 
convention that Gegenbauer polynomials of negative order are equal to zero. 



• Obviously, Eq. (2.28) is the analytic continuation from Euclidean spacetime — where the Gegenbauer poly- 
nomials obey the well-known orthonormality relations within the region — 1<2:<1 — to the Minkowski 
spacetime where the arguments of the square roots and of the Gegenbauer polynomials, depending on whether 
and/or x"^ are space-like or time-like, may be imaginary and, furthermore, could take values outside that 
region. This, however, has no influence on the validity of the above result since these polynomials being entire 
analytic functions are well defined on the whole complex plane. 



• Notice that the second term in Eq. (2.2S) is due to the exterior x— dependence of the scalar operator. Being 
proportional to x^ it vanishes on the light-cone. Furthermore, on the light-cone only the highest power z" of the 
Gegen bauer polynomial C^{z) survives and, as it should be because of ( ^.25 ), leads exactly to the expression 
( ^.24D with x^ x\ 



9 



• Finally, according to the a bove mentioned interpretation of q as some kind of 'operator symbol' within the 
Fourier integral, Eq. ( ^.2^ ) operationally has to be understood as if the polynomial inside the curly brackets 

had been written in terms of D instead of q and inserted into the operator ■!/)(O)7'''0(O). 

Now, let us remind the generating function of the Gegenbauer polynomials (see, e.g., Ref. ||2^, Eq. II. 5. 13. 1.3), 

T.^^C:{z)^T(y + i) v/T^) J.-1/2 (av/r^) e-, (2.29) 



where {2v)n = 2v{2v + 1) . . . {2v + n — 1) = T{n + 2i^)/r(2i^) is the Pochhammer symbol. Choosing z — (qx)/ ^ q^x"^ 
and a — \^ q^x"^ the local scalar operators of twist-2, Eq. ( 2.28 ), can be summed up according to Eq. (2J) to the 
bilocal scalar operator of twist-2: 

O'-^^t^x, -Kx) = ^ 0^.fi(x) U^q {^J^i^){q) [x^ (2 + xd)~l ii^q^x^} (3 + xd) 

71=0 ' 

X [^Viq^y " q^x^) ''^'j3/2 {^Vii^r-i^x^) e'"^"/'- (2.30) 

Here, the homogeneous derivations (c + xd) are required to compensate for some extra factors (c + n) which are 
necessary in order to by able to introduce the Pochhammer symbols in the denominator. Obviously, for the second 
term in Eq. ( p. 28 ), after shifting n — > n + 1 in the series over n, only one additional factor is required. — Also here q 
has to be considered as a symbol replacing the covariant derivatives sand wiched betw een the quark operators. 

Now, let us discuss the matrix elements of the twist-2 operators ( p. 28 ) and ( 2.3C ) starting with the forward case 
Pi = P2 = P- The g— integral representation of the local twist-2 operator, observing the symmetry in q and x of the 
twist-2 projection of (gx)""*"^, 

' qx 



q^x'^ 



V^:Ux',d'){xq)-+' = 

may be rewritten as follows 

= [p^^U^^d'){x^- . . J {i,r^)iq) -^[V^^l^iq^d'^){q'^^ . ..q^"+^)}. 

Obviously, the integrand is an irreducible tensor operator of order n + 1 . Its forward matrix elements are given by the 

(2) 

reduced matrix element /„ times the irreducible tensor written in terms of the momentum P-|_ = 2P. As the result 
we finally obtain 

p(2) 



{P,S\0'-^l{x)\P,S) = J^^{u{P)ru{P)) ci^. 



Px 



/•''i,..,{(v^)»cj(-|L,)-0!(v^v"-v. r - 

Now, as usual, the reduced matrix elements are considered as the moments of a distribution f'^^^z), 

/(2)= /''dzz"/(2)(z). (2.32) 



^^^-^[^n- (2.31) 



Taking this into account we obtain, after summing up over n, the following expression of the forward matrix elements 
of the non-local twist-2 operators 

{P,S\0''"\kx,-kx)\P,S) ^ u{P){x-f)u{P) J dz/(2)(z) | (2 + ^9) - Ikz-^^^ | {3 + xd) 

' (2kz^/{Px)^ - P^x^^ J3/2 (^KZy/{Px)^ - P2a;2^ ^iK.Px^ ^2.33) 



X Vtt 
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which coincides with the expression ( 2.20| ) 



More generally, when non-forward matrix elements are taken the argumentation is almost the same. But now, in 
the case of the local operators, there occur n + 1 different reduced matrix elements, /n,m, related to the monomials 
pm pn-m ^ < m < Ti which are obtained from (P2 + Pi)", and the corresponding non-forward matrix element reads: 

{P2,s2\o'z-!{^m,s^) = ('^)/(%(^2(P2)7M"(A))nfii(^',9'){^^(^Ar(a;P2)"-"}. (2.34) 

m=0 ^ ^ 

Now let us rewrite the reduced matrix elements /!% as the double moments of some double distribution f^'^\zi, Z2), 

fi% = C dzi dz2zr4""/^'n^i,^2). (2.35) 



After resumming w.r.t. n and rewriting f^'^\zi, Z2) = F^\z+, zS) we finally arrive at the expression ( 2.19| ) for the 
non-forward matrix element of the nonlocal operator O'^ '^{kx, —kx). 

Obviously, this second approach is more appropriate since, as an intermediate step, also knowledge about the local 
operators of definite twist is obtained which is related to the corresponding moments of the double distributions. 
From this approach it becomes also obvious that the mass corrections of different physical processes are related to 
each other if they can be traced back to the same operator content. 

In the following we use this polynomial technique to determine the local as well as nonlocal quark operators of 
definite twist. Thereby, we generalize Nachtmann's approach not only to non-forward matrix elements of (non)local 
operators. In addition we consider also the case of more general tensor operators having nontrivial symmetry types. 

III. TOTALLY SYMMETRIC NONLOCAL HARMONIC OPERATORS OF ANY GEOMETRIC TWIST 

The determination of the complete twist decomposition of non-local operators into the infinite tower of harmonic 
operators of definite twist is possible only for the case when their rt— th moments are related to irreducible tensor 
representations of the Lorentz group having symmetry type [n], i.e., being completely symmetric and traceless. The 
corresponding group theoretical background underlyi ng t his polynomial technique has been formulated by Bargmann 
and Todorov There, also the projection property ( |2.6| ) is proven by explicit construction. In order to demonstrate 
the efficiency of the polynomial technique and because of their physical relevance w e app ly it t o th e complete twist 



decomposition of the scalar operators N{x, —x), M{x, —x) and 0{x, —x), Eqs. (O), ( 1.12 ) and (L8), having minimal 



twist T = 3 and t = 2, respectively, as well as to the related vector and tensor operators. At the same time the 



general aspects of the twist projection, Eqs. (2.4) - (|2.6|), together with the related comments, will be exemplified. 



A. Scalar harmonic operators without external operations 

To begin with the simplest case let us study the local operators Nn (x) of degree n which are generating polynomials 
of symmetric tensors of degree n. The general formula for the decomposition of the local scalar operator Nn{x) into 
its harmonic polynomials of definite twist r = 3 -t- 2j, j = 0, 1, . . . , reads ||^ 

'"' (-iy(n + l-2j)! „,3+2i). 



where the harmonic operators of definite twist are defined by (see also Ref. [Q, Eq. 9.3.2(3)) 



they satisfy 



nNTJ^;'^\x) = Q. 



Therefore, the polynomials of Eq. (p. 2D span the space of homogeneous harmonic polynomials of degree n — 2j. 
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After substituting iV„(x) = J dq N{q){xq)" into and observing 



(n-2fc-2j)! 

we can rewrite these harmonic operators in terms of Gegenbauer polynomials as follows: 



N. 



tw(3+2i) 
n-2j 



(x) 



{n - 2j)! 



dqN{q)q^^ (^l^^ 



n-2] 



qx 



(3.3) 



Here, the factorials in front of the integral are normalization coefficients and the integration introduces a superpo- 
sition of the Fourier transforms N(q) times q^^ with coefficients being (two-sided) harmonic polynomials, 

/ii-2,(9N) = with D,hl_,^iq\x) = = D,hl,_^^{q\x). (3.4) 

The resummation of the moments Nn{x) according to N(x,—x) = '^'^^qI^^"' /nl)Nn{x) is obtained by using the 
well-known integral representation of Euler's beta function, 



B{n, m) 



r(n)r(m) 



/ dii"-i(l -t)'"-\ 

^0 



(3.5) 



Tin. + m) 

The (infinite) twist decomposition of the nonlocal scalar operator N{x, —x) into nonlocal harmonic operators 



tw(3+2j) 



jYtw(3+2j) — x) of twist r = 3 + 2j, j = 0, 1, 2, . . . , reads (remind that the summation over n for N^_2j 
at n = 2j), 



(-l)%2j 



N{x,-x)=N'-^{x,-x)+J2 i„ / 



starts 



(3.6) 



with 



7Vt-(3+2j)(x,-x) = dqNiq) {~qy {1 + qd,) (^iqx)^-q^xA "ji/2 (^^/{^xf^^^) (3.7) 



-1/2 



Let us remark that, for notational simplicity, in Eq. ( |3.6| ) we wrote the nonlocal harmonic operators for j > 1 without 
including the integration over t which is due to the normalization of the local harmonic operators, cf. Eqs. (3.1) and 
( |3.3|) . The normalization of the nonlocal harmonic operators may be read off from the terms in front of the i— integral 
in Eq. ( |3.6| ). Let us point also to the remarkable fact that the property of harmonicity is independent of j; this obtains 
immediately from the fact that the following expressions are (two-sided) harmonic functions, 

UM^) = ^ (VM' - '?'2;2) Ji/2 (I v/(9a^)' - g'a^') e'«"/' with Hi(g|x) = = Hi(<z|a:). (3.8) 

Furthermore, we observe that the factor (1 -I- qd/dq) in Eq. (B.7) may be changed into (1 -I- xd/dx) after which it can 
be taken outside the g— integration and, in Eq. (^), could be changed into (1 +td/dt); if desired, this could be used 
for a partial integration. 

The resulting expression (3.7) can be rewritten in a form which, in the case j = 0, has been already introduced by 
[Ea,R3[. Namely, using 



M-i)V!(^^)"(^) 



(3.9) 



we obtain 

^tw(3+2j)(^^_^^ 



dqN{q) {~qy 
dqN{q) {~qy 



1 



^2^,2'^ gi92;/2 



2 — q'^x^ 



qdq (^sin [^^J{qxY - q^x"^ 

qx + ^{qxf - q^x^) e(i/2) G^+VC?-)^-?^-^) 



- {qx - ^{qxf - g2x2) e(i/2)(9--V(9-)'-9'-') 
Analogous results may be obtained also in the more complicated cases which will be considered below. 



(3.10) 
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B. Scalar harmonic operators with external operations 



Now, let us consider the case of 0(x, —x) = x°'Oa{x^ — x), w here an additional power of a; occurs through contraction 
of the vector operator with a;". Therefore, the formulae (3.1) and (3.2) are to be written down for the scalar operator 
On+i(x) = J dqXfj,0^{q){xq)" , i.e., by replacing n hy n + 1 within these expressions for the decomposition into 
harmonic operators of any twist t — 2 + 2j. Then the decomposition of the scalar local operators On+i{x) reads: 



On+l{x) = 



J=0 



(-ip(n + 2-2j)! otw(2+2i) , . 



(3.11) 



where the local harmonic operators of twist r = 2 + 2j are given by 

.tw(2+2j). 



(-l)'=(n+l-2j-fc)! 



^2fen'= (□^0„+i(x)) 



4'=fc!(n + 1 - 2j)\ 



(3.12) 



After performing the differentiations □'"'^^0„+i(x) we can rewrite these harmonic operators in terms of Gegenbauer 
polynomials as follows: 



tw(2+2j), ^ 
l-2j \-^) 



(n + l-2j)! 



dgC"^(g) 

"-1-2J 



n-2i 



^2 



qx 



y/q^x 



(3.13) 



n-l-2j 



qx 



n+l-2j 



^1 



qx 



where the additional terms proportional to occur because of the appearance of derivatives with respect to (xj). 
The change in the order of the Gegenbau er polynomials and, consequently, in the order of the Bessel functions and 
the accompanyi ng fa ctors ^ (qx)"^ — q^x'^ occurring below, is also due to that fact. Of course, for j = we re-obtain 

the expression ( 2.28 ) 

Resumming with respect to rt ( > 2j ± 1 or 2j, respectively) and, again, using the representation (3.5) of Euler's 
beta function we obtain the following (infinite) twist decomposition 



0{x,^x) = 0^^'^{x,-x)+Y^ 



{-lyx^i 



dt t (1 - ty-^ o*^(2+2j)(i2.^ _^^^y 



(3.14) 



with 



Otw(2+2i)(^ fdqO'^{q)i-qyf (x^(2 + q^,) - | ig^x^) 



X (3 + qd,) (v/(<Zx)2 _ ^2^2) jg/^ [yiqx)^-q^x^) 



2j'-^il + qd,)[^iqx)^^q^X^) Ji/2 (l^i^xf^f^) 



,iqx/2 



(3.15) 



Again, by convention, the resummed nonlocal harmonic operators (tx, —tx) of twist r = 2 + 2j, j — 

0, 1,2, • • • , are obtained by changing any x — s- tx. Also here analogo us co mments are in order as for the simpler 
case of the operator iV(x, —x). Harmonicity of the nonlocal operators ( 3.15 ), □O'™*^^"'"^^-' (a;, —x) — 0, holds by con- 
struction and , again, it is independent of j; namely, it holds separately for the last, j-dependent term due to the 
equality (3^) and, therefore, mu st be fulfilled also for the combination of the first and second term. Obviously, for 
j = we recover the expression (2.30) . 

Secondly, let us consider the scalar operator M[x, —x) = x''d^M[^^]{x, —x) resulting from the skew tensor operator 
M[fj^i,]{x, —x). This operator is governed by totally symmetric tensors since the skew tensors M[^^]„ having symmetry 
type [n + 1, 1] are composed with another skew tensor x^'^d^^^ having symmetry type [1, 1] and the Clebsch-Gordan 
series of their direct product consists only of one symmetry type, [n]. 

The local operators are decomposed according to formulae (B.l) and (3.2). However, because of the special 'external' 
structure the resulting expressions simplify considerably. Namely, after performing the differentiations 



D'^+^d^a^^x^xqr = (n!/(n- 1 - 2j)!) (9^a^.a:'^)9'^'(g')'(<Z.T)"-i-2fc-2,^ 
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we arrive at the following harmonic operators in terms of Gegenbauer polynomials: 

Now, resumming with respect to n ( > 2_7 + 1) we finally get the infinite twist decomposition: 



M{x, -x) = M'-^ix, -x) + ^ ty-' M*-(3+2j")(te, -tx), (3.17) 



4^j!(j - 1)! Jo 

with 



X (2 + (3 + qd,) [^{qxY - q^x^) '^'jg/a [yiqxf-q'x^) e'^^/^. (3.18) 

Again, harmonicity of the nonlocal operators M*™^'^+^'^^ (x, — x) of twist r = 3 + 2j, j = 0, 1, 2, • • • , is fulfilled without 
taking care of j. 

C. Harmonic vector and tensor operators related to totally symmetric local operators 

The prescribed procedure may be used also for the twist decomposition of any totally symmetric (tensor) operator 
with an arbitrary number of free tensor indices which, by construction, are obtained from the corresponding scalar 
operators simply by applying appropriate derivatives with respect to x. 

The simplest of these operators, whose tower of infinite twist part starts with r = 2 and which contains the leading 
contributions to virtual Compton scattering is the operator 0^(a;, — x). Its symmetry type S is characterized by the 
Young frames [n + 1], 1 < n < oo (cf. also Ref. p| whe re we denoted that symmetry type by (i)). Its moments. 



0^„(a;), are obtained from the local scalar operators (3.11) of degree n + 1 by applying (l/(n + 1))9q. For the moment 
leaving aside the factor (n + 1)^^ we get 

[li+i] [i±i] 
d O ^Ax)- V (-l)'(" + 2-2j)! 2(j-i)otw(2+2,). . (-l)-'(n + 2-2j)! / ,^(2+2,) 

ao,Un+i{x}~ 2^ 4J(j-l)l(n + 2-j)l " ^«+i-2jW+2^ 4Jjl(n + 2-j)l \^<^'-'n+i-2j K^} 



i±ii 



Using the representation (3.5) of the beta function both sums may be rewritten as 

1 {-iy{n - 2j)\ ^tw(4+2,-). . _ 1 ^^ \-iyx^= C mil tv c'^^^+^'\tx) 



t ^^Z"^^ ^^o^^cC^ix^ ^ d.o-K.) + t ^§^,1 d..(i-t)-^(a„o-t%f (*.)), 

respectively, where 0,^T!i^2fH'^) given by Eq. (3.13) with j ^ j + 1. From this it becomes obvious that, beginning 



with twist-4, one obtains two different contributions of the same twist, namely a vector and a scalar part. 

Now, putting together these terms and resumming over n, thereby representing the normalizing coefficient l/(n+ 1) 
as Jq dAA", we obtain 



^0 



f dX 

Ol{x, -x) ^da —0{Xx, -Ax) 
Jo A 

= dA I (a„ 0'-2) (Ax, -Ax) + 4'],^^); 1^ i (1 - ty-' (9, 0'-(2+2,)) (^^^^ 

-I^"E i^l)^ X dU(l-i)-'0*-(4+2^)(Atx,-Ate)|, (3.19) 



14 



where O^^^-^+^J) (a;, -x) = 0*^(2+20+1)) (^j^ and 0^^'^'^+'^^\x, -x) are given by Eq. dslsj ). Both the integrals over 
A and t are due to normahzations and, therefore, it depends on the personal taste if they are to be included into 
the definition of the harmonic oper ators of twist r or if they are considered as part of the twist decomposition. It is 
interesting to note that in Eq. ( |3.1S| ) only the twist-2 vector and the twist-4 scalar operator survives on the light-cone. 
All other higher twist operators for j >\ are cancelled due to the factor x"^ . 

In the same manner we are able to define the twist decomposition of (completely symmetric) tensor operators which 
result from the expression 

nS - (-lp(n+l + r-2j)! 2j ^tw(2+2j) , . 

("i"--"-)"^ ^ " + (" + r) ^ 4Jj1(n + l + r-,7)! ^n+r-2j K^) ■ 

The general procedure is obvious from the consideration of the vector operator. Completely analogous to that case 
with only one derivative we would obtain a finite number of different operators of a given twist r. 

However, in the case of operators having more complicated symmetry type, e.g., being characterized by the symmetry 
type [n + 1, 1], which are necessary for the consideration, e.g., of skew symmetric tensor operators like M[c^](a;, —x), 
a theory which generalizes the formalism of Ref. is not available up to now. There exist only results for the lowest 
values of r which, however, are sufficient for the consideration of the physically relevant cases. They will be considered 
in the next chapter. 



IV. HARMONIC QUARK-ANTIQUARK OPERATORS OF GEOMETRIC TWIST 2 AND 3 

In this section we determine the harmonic quark-antiquark operators up to twist r = 2 and 3 which are relevant for 
the generalization of the parton distributions of definite (geometric) twist |^ to the double distributions appearing in 
the parametrization in the non-forward matrix elements of the corresponding non-local operators off the light-cone. 
Unfortunately, in the case of non-trivial Young frames where, up to now, no general group theoretical study being 
comparable with that of Ref. However, harmonic twist operators which are sufficient for the physically relevant 
off-cone operators have been obtained in Refs. H,^. They allow for the consideration of the vector and skew tensor 
operators up to twist 3. According to the procedure of Sect. II. C we start with their moments, i.e., the local off-cone 
harmonic operators of definite twist. These (pseudo) scalar, (axial) vector and skew tensor operators have been already 
determined (implicitly) in Ref. Appendix B. They are obtained by applying onto the harmonic tensor operators 
- which ensure tracelessness - the appropriate differential operators which take care of the required symmetry types. 

A. Local harmonic operators: Series representation 

To begin with we introduce the local harmonic operators of geometric twist-2 and twist-3 for both the chiral-even 



(axial) vector operators, 

0„„(x) = V;(0)7a(a^5)"V(0), (4.1) 

05an(x) = Vi(0)7a75(a:B)"^(0), (4.2) 
and the chiral-odd scalar and skew tensor operators, 

7V„(a:) =^(0)(a:l5)"^(0), (4.3) 

M[„^]„(x) = ^(0)a„^(a:5)"V(0), (4.4) 
together with their related scalar and vector operators 

On+i{x) = x"Oa,nix), M^n+i{x) = x^M^^^^x), M„(x) = 9''Af^„+i (a;). (4.5) 
First we consider the vector operators Oan {x) which allow for two different symmetry types [p|q| . The irreducible 



vector operators 0^^^{x) of Lorentz type (^^^^, ^^^) have symmetry types [n + 1] and can be taken over from the 
previous section. On the other hand, the twist-3 vector operators transform according to (-^^y^, ■^^-j-'-) ® ^^^) 
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and have symmetry types [n, 1]. They may be summed up to the lowest term of another infinite tower of non-local 
operators of twist t = 3, 4, . . . . The conditions of tracelessness for harmonic vector operators are as follows: 

OO'Z"{^)^0, a"OlV(x) = for T-2,3. (4.6) 

For the sake of completeness and in order to introduce a new const ructiv e element let us reconsider the twist-2 
vector operator. It is obtained from the scalar twist-2 operator, cf. Eq. (3.11), 

r)tw2 ( _ (-1) {n+1- fc)! 2fe pfc^ / X 

k=0 ■ ^ ^' 

by means of a partial differentiation, da, together with the normalization factor l/(n + 1) as follows: 



— T 

+ 1 ^ 

fe=0 



(Ji|^.-P.(fc)tfO„+,(x). (4.7) 



Here, we introduced the operation 



Va{k) ^ {k + 1 + xd)da - I XaD (4.8) 



which is a generalization of the inner derivative (for the latter, cf. Ref. [|| and, more detailed, Refs. ||,0) off the 
light-cone and its extension to arbitrary values of fc > 0. (A similar derivative in momentum space already has been 
introduced in Ref. |p^ .) As a side remark we mention that with the help of this generalized derivative also symmetric 
tensor operators of rank 2 and higher can easily be construct ed. 

The twist-3 vector operator which satisfies the conditions (4.6) is, with the convention ^[^i/j = ^ (^^ly — 



k=0 



The twist-4 vector operators having symmetry type \n + 1] can be read off from the general result of Sec. III.C. 
However, there occur also twist-4 vector operators having symmetry type [n, 1] which have been determined on the 
light-cone only H,^. 

The scalar twist-3 operator is of Lorentz type (^, ^) and also has been considered in Sec. HI; it reads 

^'^^(^)-E ^"Tw'^' ^"°'^"(^)- (4.10) 

fc=0 

Now, let us give the decomposition of the skew tensor operator M^ai3]n{x) into its twist-2 and twist-3 parts. This 
determination is simplified by the observation that these operators are related to the vector operator M^n+i{x) — 
x^M^^^^nix) and the scalar twist-3 operator M„(x) = 9^M^„_|_i(a;), respectively. The twist-2 tensor operator 
■^[ap]ni^) transforms according to (^^-j^,f) © (fj^^^); whereas the twist-3 operator which is obtained from the 
trace terms of Af*^j^(x) is of Lorentz type (2-, ^). They have symmetry types [n + 1, 1] and [n], respectively, and are 
given by 



2 
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and 



with 



(n + 2) 



(4.12) 



n 



k=0 



k)\ 



1 

n 



k=0 



(-l)'=(n-l-A:)! 
4'=fc!n! 



a;2'=I?/3(fc)n''M„(a;) 



(4.13) 



Besides the twist-3 operator (4.12) resulting from the trace terms of the twist-2 operator (4.11) there exists also a 
genuine twist-3 operator having Lorcntz type ^) and being governed by the symmetry type [n, 1, 1]. It is given 
by the following expression (with the convention Aap-y = ^{A^^p^] + Ap^^^] + A^^^p-^)): 



k=0 



n + 2 



E 

k=0 



x^'^x-^{s^JfjV,^{k) ~-5^^xpd^p-'ik)}n^M^^,]„ix). (4.14) 



By construction, Eqs. (4.11), (4.13) and (4.14) are (traceless) harmonic polynomials obeying 
□A^M]n(^) = 0, d"Ml:p]^{x) ^ = dPMl:p]^{x), T = 2, 3, 

and 



a 71 ,.rtw3 



(4.15) 



(4.16) 



Let us point to the structural similarities of the operators (4.7) and (4.13) on the one hand and ( |4.9| ) and (4.11) 
on the other hand. - Furthermore, we observe that when organizing the various expressions such that everywhere 
the partial derivatives appear to the right of the fcth power of x'^ the differential operations ensuring the symmetry 
type are transformed according to da — > 'Da{k) whereas the generalized (first order) inner derivative x^^'^p^ remains 
unchanged because it commutes with x^ . This might be used as a hint how to extend harmonically the operators 
when they are given on the lig ht- cone in terms of in ner de rivati ves d^ = 2^/3(0) and ^[qQ/j]. Namely, a short look 
on our results, Eqs. (4/7), ([4.6|), ( 4.10 ), ( 4.11 ), ( 4.12 ) and ( 4.14 ), shows that the expressions on-cone are obtained 
by restricting the sums to those terms with k — Q and replacing x — > i, da ^ d/dx"^. Let us mention that this is 
nothing else but another form of the harmonic extension which is well-known in group theory HJI^,^ . 

In the case of symmetric traceless tensors the harmonic extension is uniquely defined and proved in Ref Q . In all 
the other cases it can be used heuristically according to the above observation. This harmonic extension gives an 
one-to-one relation between homogeneous polynomials on the light-cone and the corresponding harmonic polynomials 
off-cone. It is important to note that the unique harmonic extension must not destroy the type of the Lorentz 
representation. 



According to Eqs. (4.7), (4.9), (4.10), (4.11), (4.12) and (4.14) we may introduce the orthogonal projection operators 
onto the subspaces with definite spin and, therefore, with definite geometric twist r = 2, 3 as follows: 



Oir) 



'S)L(^)=^ir^0(5),„(a:), 
N^:\x)=V(:^Nr.{x), 



(4.17) 
(4.18) 
(4.19) 



The corresponding twist projectors onto the light-cone are discussed in more detail in Ref. 
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B. Resummed harmonic operators with twist r = 2, 3 



In this subsection we rewrite the harmonic operators of definite twist t — 2,3 which have been given in Subsect. A 
in terms of Gegenbauer polynomials analogous to the presentation in Sect. II for the scalar harmonic operators. In 
a second step we resum the infinite series (for n) of Gegenbauer polynomials to the related Bessel functions thereby 
obtaining nonlocal harmonic operators of (the same) definite twist. 

For notational simplicity we introduce the following abbreviations for the homogeneous polynomials and their 
related functions: 



qx 



(4.20) 
(4.21) 



Using the following functional relations of the Gegenbauer and Bessel functions [ pTj , 

^(^-Va(z)) --^-Va+i(z), 
2^Ja(z) = -\J\[z) + zJa-i(2;), 
respectively, one obtains the following results for the partial derivations of these functions: 

daKXq\x) = V [qaK^_\{q\x) - x oKn-2{(l\x)^ , 
da'Hu{q\x) = [\qa{'^v +l + xd)- \{iqfxo^ Hi,+i{q\x). 



(4.22) 
(4.23) 



The local twist-2 vector oper ator is obtained through derivation with respect to of the local twist-2 scalar 
operator as follows, cf., Eq. (4.7), 

or (2^) = (^^^ J dS {^l,i^){q) {x^ hl{q\x) ~ \ q^x' hl^M^)] 



A^q {i^l^.^){q)iy5^^ hl{q\x) ~ q^^x^ hl^M^) (4-24) 
+ 2a;^g„ hl_^{q\x) - (x^'x^g^ + q^q^x^) hl^^(q\x) + \ q'^Xa.x^q^ hl_^{q\x) 



and the corresponding resummed nonlocal twist-2 vector operator reads 

0':'^X,-x)=d„ f dt fd\{{pj^ij){q){x>^{2 + qdg)^litq^x^}{3 + qdg)n2{q\tx) 



d''q{i,^,,i,){q){2 + qdq) dt\ {'i + qdq) 6^ - itq^x^ n2{q\tx) 



(4.25) 



[i + qdq) ({A + qdq)itq^x^ - \ (itf {q^ X^' x ^ + x^ q^ q^)^ + \{ltf q^" q^ x o,x''\nMix) 

Here, the equality l/{n + 1) = jj^ dt has been used; analogous equalities will be used in the following. 
The resummed local twist-3 vector operator reads, cf., Eq. (4.9): 



{n + lf 



d^q{i,-1^i:){q)xP}^2{n+l)6lqp^ hl_,{q\x) + {n + 2) x^J^^q^ hl_^{q\x) (4.26) 
- 4x[„(j0]g'' hl_2{q\x) + 2x[aqii]x''q^ hl_^{q\x) 
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and the bilocal twist-3 vector operator is given by 

= / d^g(^7^V)(g) / dt {2 + qdg) x4 [2 iUf„g^](2 + gc),) - {it)H^^xp]q']n2{q\tx) (4.27) 



itX[aqi3] 



2Ytq^'{b + qdq)-{Ytfx^'q^ Hsiqltx) 



Let us remark that in the expression (4.27) by partial integration the term containing (1 + tdt) contributes only 
throu gh it s surf ace te rms for i = 1. Remind also that according to the special structure of H^{q\tx) = Hu{tq\x) in 
Eqs. (4.25) and (4.27) the homogeneous derivation qdq may be replaced by either xdx or tdt- Analogous observations 
could be made below; however, there the a;— derivatives in the local as well as the nonlocal case will not be taken 
explicitly. 



The local twist-3 chiral-odd scalar operator will be given only for completeness, cf., Eq. (3.3) 

K^^x)^ J d4g(^V)(9) Kiq\x), 
and for the corresponding bilocal one, one gets 

N''"^ix,-x)^ f d^q{^P4'){q)ii + qd,)n,iq\x). 



(4.28) 



(4.29) 



The local resummed twist-2 chiral-odd skew tensor operator is given by, cf., Eq. (4.11), 
^[at]n(^) = ;^{'53]-^^[o5/3]5^} / d\{^a,^^){q){x'^ hl{q\x) - ^q-^ x' hl_,{q\x)}, (4.30) 
and the related bilocal twist-2 skew tensor operator reads 

Mf„-|](a:, -^) - 2y / dt t (1 + qdg){{2 + qdq) 6^^dp^ - x^^dp^d^^] [x" (3 + qd,) - \itq- x^'\n2{q\tx) . 

(4.31) 

For the related, resummed local twist-3 skew tensor operator resulting from the trace terms we obtain 

2 



^ ^' n[n + I) 



'^{adfj\ \ i'^q {■'p(^f^uip)iq) q'^x" hl_-i^{q\x), 



(4.32) 



leading to the following twist-3 bilocal skew tensor 



K"l](^'-^) = 2^["^/3] d\{^a^,i;)iq)iq^x''{2 + qdq) / dtH^iqltx). 



In addition, the independent twist-3 skew tensor having symmetry type [n, 1, 1] has the following local form: 



(4.33) 



M[:^]ni^)^^x'\S^JP^]--5i^^xpd,^d-^\ I dS{^a,,^)iq)K{q\x) 



leading to the following expression for the twist-3 bilocal skew tensor, 



(4.34) 



Ml:^^ix,-x)^Sx''\d^JP^]xd-2d^^X0d^]d''\ I dS{^a^,,4'){q)il + qdg) j dt^^ HMtx) . (4.35) 
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V. POWER CORRECTIONS OF NON-FORWARD MATRIX ELEMENTS 



A. Double distributions from quaric— antiquark operators with definite twist r 



According to the different approaches exemphfied in Subsect. II. B and II. C in order to obtain the non-forward 
matrix elements for the quark-antiquark operators with definite twist r we have only to take matrix elements of the 
expressions given in the last Subsect. IV. B and to use formula (2.11) of Subsect. II. A. Let us emphasize that the spin- 
dependence of the hadronic states \Pi,Si),i = 1,2, is completely contained in the independent kinematical structures 
/Cp(P) being relevant for the processes under consideration, e.g., the Dirac and the Pauli structure u{P2, 82)1 uu(Pit Si) 
and u(i-2, S2)(T^yP':Lu{Pi, Si)/M, respectively, in case of the virtual Compton scattering (see, e.g., p5|). 

Let us now give the final expressions for matrix elements of those nonlocal operators which are relevant for the 
power corrections of the various processes. The corresponding expressions for the local operators are easily obtained 
by taking moments. 



For the non-forward m atrix element of the nonlocal twist-2 and twist-3 vector operator, Eqs. (4.25) and (4.27) 
obeying the replacement ( ^.11 ), one gets 



{P2,S2\0'J^{kx,-kx)\Pi,Si) ^IC^iP) J BZFi'^\Z)daJ dt |a;^ (2 + xS) - ^ i Kt (P^Z) | (3 + xS) H2(PZ|Kte) 

= /C;i(P) J DZFi2)(Z)^ di(2 + .Ta) I {3 + xd)d^-iKt{P''Z)xa H2(PZ|Kte) 
{3 + xd) (^{4 + xd)iKt (P„Z)a;'' - liint)^ {{PZ^x^Xa + (P^Z)(P„Z) x^)) 

+ jiintf {¥''Z)(VZfxo,x^]n3{VZ\Ktx) 



(5.1) 



and 



{P2,S2\0':'^Kx,-Kx)\PuSi) = 2/C°(P) J DZFf ^ ^'^{^[Jp] {l + xd) ~ X[^dp]d^}ni{PZ\Ktx) 

= /C;i(P) / DZF^iZ) {2 + xd) I dtxl^l 2iKt5'^^{Pp)Z){2 + xd)~{iKtf5'^^xp]{PZf 7^2(PZ|Kte) 



- iKta;[„(P^]Z) 2iKi (P^Z) (5 + xa) - (iKt)2x^(PZ)2 7i3(PZ|Kte)|. (5.2) 

In the case of the axial vector operators 0'^^{k.x^ —kx) the kinematical structures /C5^(P) contain an additional 75 

and the double distributions should be denoted by G'"J\z) but otherwise the structure will be e xactly the same. 
The non-forward matrix elements of the nonlocal twist-3 chiral-odd scalar operator, Eq. (4.29), reads: 



(P2,S2\N*'"^{kx,-kx)\Pi,Si) =/C"(P) f T)ZE\^\Z) {I + xd) ni{PZ\Kx). 



(5.3) 



The non- forw ard m atrix elements nonlocal twist-2 chiral-odd skew tensor operator and of the related twist-3 oper- 
ator, Eqs. (4.31) and (4.33), are given by 



(P2, 52|M[!r|]('*a:, -Acx)|Pi, 5i) = 2/Cf^,](P) j DZH^iZ) j^^tt{l + xd)[{2 + xd)5ldp^ - x^^df^^d^^] 

X x" (i^ xd) - \ int {rZ)x^ 'H2{PZ\Ktx), 

(P2, ^2|M[a|]('«a^- -Kx)\Pi,Si) = 2/Cf^,](P) J DZi/P(Z) X[^dp] m (P^Z)a;''(2 + xd) dt'H2{n\Ktx) . 



(5.4) 



(5.5) 
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Finally, the non-forward matrix element of the independent nonlocal twist-3 skew tensor, Eq. (4.35), reads: 



P2,S2\Ml:;^]i^x,-Kx)\P,,S,) = 3/Cf^,](P) / BZHPiZ)x-y[s^J-^f,d^^xd-2 6i^^xpd^^d''} 



x{l + xd) / dt niiPZ\Ktx). 

Jo 2t 



(5.6) 



The matrix elements of the scalar and vector operators are obtained by multiplying the above expressions by x" (of 
course, for the operator M^^^^{kx, —kx) no vector operator exists because of its symmetry type). For completeness let 
us als o give the p ower corrections of the double distributions of higher twist in the scalar cases. Using the expressions 
(131), Q and (131) we obtain 



{P2,S2\0''"^^+^^\kx,~kx)\Pi,Si) =JCliF) / DZ F^^+sj) (Z) (iPZ) 



Xf" (2 + xd) - i iK (P^Z) X^ (3 + xd) 'H2{n\Kx) - 2j^^^ (1 + xd) 7il(PZ|K.T) 



(P2,5'2|iV*^(='+2^'^(ra,-ra)|Pi,S'i) =/C''(P) / DZ£:f+2^)(Z)(iPZ)2J (l + a;a)Hi(PZ|Ka;) 



(5.7) 



(5.8) 



and 



5'2|M*^(3+2j)(Kx,-Ka;)|Pi,S'i) = /Cf^^](P) / T>Z H'^J^+'^^^I) lV^"Lx'' {in)^^ [2 + xd) {2, + xd) 7^2(PZ|kx), (5.9) 



respectively. Let us remind that, eventually, the ^-integration which appeared, e.g., in Eq. ( |3.6| ) should be included 
into the definition of the operators of definite twist and, therefore, it could appear also here. 

Now, a couple of remarks are in order. First, as indicated, different tensorial structures of the operators lead to 
different kinematical structures. In the forward case they simplify or eventually disappear, e.g., the Dirac structures 
are to be replaced by and 5*^ for the vector and axial vector case, respectively, whereas the Pauli structures 
vanish. Second, going on-conc the double distributions remain the same. Therefore, their evolution is determined 
by the anomalous dimensions resulting from the renormalization group equation of the corresponding light-ray op- 
erators. Third, taking the forward limit on the light-cone the double distributions are projected onto the quark 
distribution functions determined in Ref. In this connection we emphasize that these distribution functions are 
uniquely determined by the geometric twist which differs from the dynamic twist being used for phenomenological 
considerations. 

Furthermore, we have to mention that - contrary to the case when the non-forward matrix elements are restricted 
to the light-cone where the decomposition of the non-local quark-antiquark operators into operators of definite twist 
terminates at finite values Tmax ~ the off-cone decomposition results in an infinite series of any twist. Therefore, the 
decomposition of the double distributions related to the undecomposed operators, /^(Z, P^Pj, a;^, ^^), results also in 
an infinite series. This becomes obvious especially in the scalar cases considered in Sec. III. 

The double distributions occur in the scattering amplitudes of various physical processes, like virtual Compton 
scattering in the generalized Bjorken region. There, after multiplication with appropriate coefficient functions Cr(x) 
a Fourier transformation with respect to x has to performed relating x to the (inverse) momentum transfer l/Q of 
the process. Thereby, the variables (xPZ) and a;2(PZ)2 are changed into (QPZ)/g2 and (PZ)V(9^- The Fourier 
transformation will be simplified when the Poisson integral is used for the Bessel functions. This will be considered 
in more detail in another paper. 

Another application of our general procedure is the consideration of the meson distribution functions which does 
not require such an additional Fourier transformation. 



B. Mass corrections of vector meson distributions 



Based on bilocal light-cone operators of definite twist-2, 3 and 4 the distribution amplitudes (DAs) of the p-meson 
have been already studied in Refs. |^,|ri|. The relationship between the DAs of geometric and dynamical twist has 
also been discussed (see, also, Ref. )7^These considerations can be extended to the other vector mesons by trivial 
substitutions. It is an advantage of this approach that operators of well-defined twist are used and, therefore, it is not 
necessary to use operator relations, like dynamical Wandzura-Wilczek relations, in order to isolate contributions of 
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twist-3 and 4. The difFerence, in the non-forward case, between geometric and dynamical Wandzura-Wilczek relations 
was pointed out in Refs. |pl| , M - 

In this Subsection we extend that analysis of the p-meson DAs to include all meson momentum resp. mass terms of 
harmonic twist-2 and twist-3 off-cone operators in order to get the power corrections to the results obtained earlier. 
Vector meson mass corrections of order were already discussed by Ball and Braun and resummed meson mass 
corrections in the scalar case was given by Ball |3^ (see also Ref. P9|]). 

The structure of the mass corrections of p-meson DAs of geometric twist are, from the group theoretical point of 
view, similar to the target mass corrections in deep inelastic scattering, which can be resummed using Nachtmann's 
method [Q. Also here, the basic tool, in order to obtain the mass corrections, are the harmonic operators of 
definite geometric twist which have been determined in Sec. IV. Obviously, they are the harmonic extensions of the 
corresponding light-cone operators which already are used in Refs. for the classification of the corresponding 

meson light cone DAs with respect to geometric twist. 

Let us now introduce the distribution functions for the harmonic operators of geometric twist sandwiched between 
the vacuum and the meson state, (0|C> p"''( x, —x')\p{P, A)). Thereby, we adopt the definitions of Chernyak and Zhit- 
nitsky |36| in the terminology of Ref. j^7| . As usual, these matrix elements are related to the momentum Pa and 
polarization vector ei^"* of the meson with helicity A, P^ = m^, e^'*'-' • e*^^-' = —1, P ■ e^'^^ = 0, m denotes the meson 
mass. 

First, we consider the chiral-even vector operator. Using the twist projections we introduce the moments of the 
meson DAs, of twist r according to 

(0|Oi;Hx)|p(P, A)) = f.mViV/ (e^'^(P2:)"<l>M) , (5.10) 

where fp is the vector meson decay constant. The corresponding meson distribution amplitudes $'•^■'(0 9.re given by 
inverting the moment integral, 



(5.11) 



These distribution amplitudes are dimensionsless functions of ^ and describe the probability amplitudes to find the 
p-meson in a state with minimal number of constituents (quark and antiquark) which carry the momentum fractions 

u (quark) and {1 — u) (antiquark). 

Taking into account expression (4.24) we obtain the local twist-2 matrix element as follows: 

(0|Or„'(x)|p(P,A)) = ^-^/p™$(2)|e(^)/i^(P|2;) + 2P„(eWx)/.Li(P|x)-m2^^ (5.12) 

which is analogous to Wandzura's expression [Q for the target mass corrections in deep inelastic scattering in the 
x-space. After resummation we get the bilocal matrix element of geometric twist-2: 

(0|Or (x, - x)|p(P, A)) = fpmj^' dt d^^(^'\0 (3 + m (2 + ^ 

X |ei^)H2(P?|te)-|- (^(4 + ea^)(i^i)P„-i(iet)Wx„^(e(^)a;)H3(P?|te)|. (5.13) 
An analogous calculation using expression ( 1.26| ) gives the local twist-3 matrix element 

{0\O'ZH^MP, A)) = ^^^^ /pm $(f)x''{2(n + 1) e{^^Pp^hl,,{P\x) + (n + 2) m\J^]hl_,{P\x) 

+ 2m^X[aPp]ie^^^x)hl_^{P\x)}; (5.14) 
the resummed bilocal matrix element of twist-3 reads 

r-l 



\0'-Hx,-xMP,X))^fprn J ^d^^^'HO {2 + ^d^) x^ }^2{iOe\yp]n2{P^\x) 

+ m2 f dt((iCi)2a;[„e^]^H2(PC|te) + (iei)3a:[c.P^](e(^)a;)H3(P^|te))|. (5.15) 
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It is well-known that the twist-3 operator O^J'^{x, —x) is related to other twist-3 operators containing total derivatives 
and oper ators of Sh uryak -Vainshtein type by means of QCD equations of motion [^6l|35| ] . Therefore, the matrix 
elements (5.14) and (5.15) include the contribution of the twist-3 operator containing total derivatives which is as 
large as those from twist-2 (see also |34j). In addition the next higher twist contributions of the chiral-even vector 
operator are of twist-4. 

Now we consider the chiral-even axial vector operator O^anix) with the same twist projectors as for the chiral-even 



vector operator Oan (x) . We define the corresponding moments of the meson DAs E^n"^ by 



(5.16) 



where fj denotes the tensor decay constant. First, we observe that the twist-2 contribution vanishes. The nontrivial 
local vacuum-to-meson matrix elements of this axial vector operator are of twist-3: 



{o\ori{xMP,\)) = l{fp~.fj 



m„ -I- ma 



2 V ■ - m 
and the bilocal matrix element of twist-3 reads 



-^13 ^ 



PfiXi, hj^[P\x), 



(5.17) 



{0\OrJ'{x,-xMP,\)) = ^ (/p - /J ^ m ej^'^ef^ P,x, dg S^^) (^ (1 + ) H, (Pg|x) . (5.18) 

By the way the twist-4 matrix eleme nt als o van ishes and the next higher twist contribution would be of twist-5. 
Also here, the matrix elements (5.17) and (5.18) include the contribution of the twist-3 operator containing total 
derivatives. 

The matrix element of the chiral-odd scalar operator is defined as 



(0|iVr(:r)|p(F,A))=-i(/J-/p 



13)((e(^)x)(Px)"Tf 



(5.19) 



/o^ 

where T„ is the moment of a spin-independent twist-3 distribution function. Using expression (4.28) the local matrix 

(5.20) 



element is given as 

m^rHxMP, A)) = -i(fj - f^!Il^±J^) (e(')^)m' Tl^^hi{P\x), 



and for the bilocal matrix element of twist-3 we obtain 



7v*-^(x,-x)|p(p,A)) = -i 



mu + md \ I (A) 



d£,T^^\0{^ + m'ni{Pi\x). 



(5.21) 



The next higher twist contributions of the scalar chiral-odd operator are of order twist-5. 

Now, we consider the matrix elements of the chiral-odd skew tensor operators. The corresponding moments of the 
wave function are introduced by 



\M[:U-)\p{P, A)) = i/JPg: ((e(^)p. - e(^)p,) (P.)"*!^ 



(5.22) 



The local matrix element of the skew tensor operator of twist-2, performing the differentiations in the expression 
(4.30), is given by 

(0|Mf-|„(a;)|p(P,A)) = i/J*(2)|_l_e[^^)p^];j2(p|^) + ^llL. ^.^J^) hl_^{P\x) 



m 



(n-K2)2(n-H 1) 



-[^x^J^:^hl_^{P\x) ~ {2x^J^^{Px)+Ax^^Pp^{e^^^x))hl_^{P\x) 



+ 12x[„P^](e(^)x)/i^_2(P|x) - &x^^Pp^{e^''^x){Px)hi_^{P\x)]y 



(5.23) 
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and the resummed bilocal matrix element of twist-2 reads 



(0|Mf-|](x, A)) = i /J j ^ d^ (3 + e^e) |2 (2 + e^c) e}^)p^]H2(Pe|a;) 



dti 



+ (2 {A + ^d^)mx[^e^^^ - (ie<)'(x[„e^-;^(Px) +2x[„P^](e(^)x)))H3(F^|te) 

+ (2 (5 + ^5^) + (iei)'a:[„P0](e(^):E) - (4 + ^9^) H4(P^|te) 



(5.24) 



The local matrix element of the skew tensor operator of twist-3 is obtained from (4.32) as follows: 

2 



(0|M[r|]„(2;)|p(P,A)> 



{n + 2)n 

and the corresponding bilocal matrix element reads 



ifJ^ifW[x[J^hl^,iP\x) + 2xi^P^^ie(^'>x)hl^,{P\x)}, (5.25) 



(0|Mf-|](x,-2;)|p(P,A)) =-2i/>2£^ f^dC¥'HOa + (d^) 

X {mx[^e^^^H2iP^\tx) + (3 + ^9?) m^X[^Pp]{e^^^x)H3{P^\tx)} . 



(5.26) 



The nex t hig her twist contributions of the skew tensor operator would be of twist-4. Again, the matrix elements ( p. 20 ), 
( 5.21 ), ( |5.25| ) and ( 5.26 ) include the contribution of the twis t-3 o perator containing total derivatives. Let us note, 
that the contribution of the additional twist-3 operator, Eq. (4.34), vanishes. 

Finally, we remark that after projection onto the light-cone the matrix elements, which have been introduced in 
Rcf. jl^, are recovered. This can be easily checked by using the expansion of the Gegenbauer polynomials or the 
Poisson integral representation of the Bessel functions. Let us also point to the essential fact that the x^— dependence 
is contained in the polynomials h'^{P\x) and the functions Ti.^{P\x) in the case of the moments and the wave functions, 
respectively. Since the latter is given by 



2r(z.)H.(P|x) = 4"^ di(l - t2).-igi(P.)(i+Vi-™^-V(P-)^)/2 



(5.27) 



one would obtain a fairly simple expansion in powers of x^ / (Px)'^ . As has been pointed out for general distribution 
amplitudes in Sec. II the meson DAs of definite twist are independent from any coordinates or momenta and therefore, 
coincide on-cone and off-cone. The whole information about the power corrections is contained in twist projection of 
the operators from which the matrix elements are taken. 



VI. CONCLUSIONS AND COMMENTS 



In this paper, we introduced a general procedure of parametrizing non-forward matrix elements of non-local QCD 
operators Or{nx, —kx) by multi- variable distribution amplitudes fa^\Z) of well-defined geometric twist r, namely, 
single variable hadron distributions, double distributions, triple distributions etc., times position-dependent coefficient 
functions c^'^^ as well as kinematical factors JCr related to the Lorentz structure of the QCD operators and the hadron 

(r) 

states sandwiching the operators Op . The procedure relies on the unique twist decomposition of non-local operators 
off the light-cone leading to an infinite series of operators with growing twist. This decomposition is completely of group 
theoretical origin and is equivalent to the decomposition of the local operators into irreducible tensor representations 
of the Lorentz group. In the case of operators whose moments are related to totally symmetric tensors, using the 
procedure of harmonic extension |^, we were able to determine the decomposition completely. For operators with 
non-trivial symmetry type, due to the lack of a general theoretical framework, only the terms up to twist t = 3 have 
been determined. However, this covers already almost all the phenomenological interesting cases of quark-antiquark 
operators. 

Using these results we determined the off-cone power corrections to various double distributions and the vector 
meson wave functions being the inputs of the corresponding scattering amplitudes and hadronic form factors, respec- 
tively. These power corrections, in the case of local operators, are expressed by a finite sum in terms of Gegenbauer 
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polynomials being multiplied with the moments of the distribution amplitudes and, in the case of non-local operators, 
by a finite sum in terms of Bessel functions now multiplied with the distribution amplitudes directly. As a remarkable 
fact, which may be read off from Sec. V, we observe that (the moments of) the distribution amplitudes are independent 
of whether the operators are taken off-cone or on-cone. Furthermore, the off-cone expressions are obtained from the 
on-cone ones by harmonic extension. 

This very encouraging behaviour strongly relies on the definition of the distribution amplitudes through matrix 
elements of operators with definite geometric twist - contrary to their possible definition with respect to 'dynamical' 
twist. This, of course, would extend a procedure to the non-forward case which has been introduced for the forward 
matrix elements (on the light-cone) in Ref. |^ . As a consequence it is not necessary to use any further dynamical input, 
like equations of motion, and to take into account operators containing total derivatives in order to get expressions 
of well-defined geometric twist: their contribution is already contained in the expressions of given (geometric) twist. 

Concerning the computation of the scattering amplitudes resp. form factors of physical relevance some Fourier 
transformation has to be carried out whose result mainly depends on the (singular) coefficient functions as well as 
on the various Bessel functions. However, modulo specific tensorial structures, one simply has formally to substitute 
X iq/Q"^ leading to the replacements of (xPZ) {qn)/Q'^ and x^{nf /{xflf {[fZf / Q'^) {Q"^ / {qfTiY , 
especially in the expressions h'^{J'X\x) and 7i^(PZ|a;). The currently most interesting example is the (deeply) virtual 
Compton scattering whose consideration has been postponed to another paper. 

There is a further point deserving attention. Since the distribution amplitudes of definite twist are given already on- 
cone their renormalization group induced Q^— evolution, which is determined by the (non-local) anomalous dimension 
of the non-local light-ray operators of definite twist, should be clearly separated from the kinematical power- like target 
momentum resp. mass corrections. 
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